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Abstract We study collective wind configurations pro- 
duced by a number of massive stars, and obtain den- 
sities and expansion velocities of the stellar wind gas 
that is to be target, in this model, of hadronic interac- 
tions. We study the expected 7-ray emission from these 
regions, considering in an approximate way the effect 
of cosmic ray modulation. We compute secondary par- 
ticle production (electrons from knock-on interactions 
and electrons and positrons from charged pion decay), 
and solve the loss equation with ionization, synchrotron, 
bremsstrahlung, inverse Compton, and expansion losses. 
We provide examples where configurations can produce 
sources for GLAST satellite, and the MAGIC, HESS, or 
VERITAS telescopes in non-uniform ways, i.e., with or 
without the corresponding counterparts. We show that 
in all cases we studied no EGRET source is expected. 
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1 Single and collective stellar winds 

LBL, WR, O and B stars lose a significant fraction of 
their mass in stellar winds with terminal velocities that 
can easily reach Voo ~ 10'^ km s^^. With mass loss rates 
as high as = (IQ-^ - lO^^) Mq yr"!, the den- 
sity at the base of the wind can reach 10~^^ g cm~'^ 
(e.g., [I]). Such winds are permeated by significant mag- 
netic fields, and provide a matter field dense enough as 
to produce hadronic 7-ray emission when pervaded by 
relativistic particles. A typical wind configuration ([1, 
[2ll3]) contains an inner region in free expansion (zone 
I) and a much larger hot compressed wind (zone II). 
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These are finally surrounded by a thin layer of dense 
swept- up gas (zone III); the final interface with the in- 
terstellar medium (ISM). The innermost region size is 
fixed by requiring that at the end of the free expan- 
sion phase (about 100 years after the wind turns on) 
the swept-up material is comparable to the mass in the 
driven wave from the wind, which happens at a radius 
^wind = T4o(3M^,/47rpoK^)^/^, where po « mpUo is the 
ISM mass density, with nip the mass of the proton and 
no the ISM number density. After hundreds of thousands 
of years, the wind produces a bubble with a radius of 
the order of tens of parsecs, with a density lower (ex- 
cept that in zone I) than in the ISM. The matter in the 
inner region is described through the continuity equa- 
tion: — Anr^ p{r)V (r) , where p{r) is the density of 
the wind and V{r) = \4o(l — Ro/r)^ is its velocity. 
Voo is the terminal wind velocity, and the parameter p 
is ~ 1 for massive stars ([1]). Rq is given in terms of 
the wind velocity close to the star, Vq ^ 10~'^Vtx), as 
Ro = i?*(l — (Vo/Vooy^^). Hence, the particle density is 
n{r) = Af*(l - Ro/ryf^ /{ATTiripVoor'^)- In the case of a 
collection of winds, a hydrodynamical model taking into 
account the composition of different single stellar winds 
needs to be adopted. As we see below, differences with 
the single stellar wind case are notable. 

Here we adopt a similar modelling to that of Canto et 
al. (2000) [3] (see also [illl[7]) to describe the wind of a 
cluster (or a sub-cluster) of stars. Consider that there are 
iV stars in close proximity, uniformly distributed within 
a radius i?c, with a number density 

Each star has its own mass- loss rate (Mi) and (terminal) 
wind velocity (Vi), and average values can be defined as 
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All stellar winds are assumed to mix with the surround- 
ing ISM and with each other, filling the intra-cluster 
volume with a hot, shocked, collective stellar wind. A 
stationary flow in which mass and energy provided by 
stellar winds escape through the outer boundary of the 
cluster is established. For an arbitrary distance R from 
the center of the cluster, mass and energy conservation 
imply that 
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where p and V are the mean density and velocity of 
the cluster wind flow at position R and h is its specific 
enthalpy (sum of internal energy plus the pressure times 
the volume), 
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with P being the mean pressure of the wind and 7 being 
the adiabatic index (hereafter 7 = 5/3 to fix numerical 
values). From the mass conservation equation we obtain 
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and the ratio of the two conservation equations imply 
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The equation of motion of the flow is 
^dV _ dP 
dR ~ ^dR 

which, introducing the adiabatic sound speed c, 

2 P 
c = 7— , 

P 

can be written as 
.dV _ 1 d{pc^) 
Zr ~ ~7 dR 
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From the definition of enthalpy and Equation ([S|), the 
adiabatic sound speed can be expressed as 
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Using Equation (O, its derivative dp and Equation ([12 
in pT|) one obtains 



dR 
~R 



dy_ 

V 



(7 - i)vJ - (7 + i)v^ 
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which can be integrated and expressed in more conve- 
nient dimensionless variables {v = V/Vw and r = R/Rc) 
as follows 



1 



57-1-1 
\ 

7-1 



-, -(37+l)/(57+l) 



(14) 



with A an integration constant. When R > Rc, i.e., out- 
side the cluster, by definition n — and the mass con- 
servation equation is 



An 



-Rc^nA'U = inR^pV, 
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where the middle equality gives account of the contribu- 
tion of all stars in the association, and Afassoc — X^i 
is the mass- loss rate at the outer boundary Rc- Substi- 
tuting Equation (fT^ and (fT5|) into the n = realization 
of Equation (fTT|) one obtains 
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and integrating, the velocity in this outside region is im- 
plicitly defined from 



w(l - w2)i/(7-i) = Br 



(17) 



with B an integration constant. Having constants A and 
B in Equations (fH)) and (fT7|) . see below, the velocity 
at any distance from the association center can be deter- 
mined by numerically solving its implicit definitions, and 
hence the density is also determined, through Equation 

or msi). 

From Equation ()17|) . two asymptotic branches can 
be found. When r — > 00, either w — s- (asymptotically 
subsonic flow) or w — > 1 (asymptotically supersonic flow) 
are possible solutions. The first one (subsonic) produces 
the following limits for the density, the sound speed and 
the pressure 
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If Poo is the ISM pressure far from the association, the 
constant B can be obtained as 
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The velocity of the flow at the outer radius r — 1 
follows from Equation (fT7|l 
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and continuity implies that 
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Equation p4p implicitly contains the dependence of v 
with r in the inner region of the collective wind. Its left 
hand side is an ever increasing function. Thus, for the 
equality to be fulfilled for all values of radius (0< r <1), 
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Table 1 Examples of configurations of collective stellar 
winds. The mass is that contained within 10 Rc, and is shown 
as an example, no is the central density. 



Model 




K, 


Rc 


no 


Wind mass 




[Mo yr-^1 


[km s"^] 


pc 


cm"^ 


[M0 ] 


A 




800 


0.1 


210.0 


0.13 


B 


10"" 


800 


0.3 


23.3 


0.39 


C 


5 X 10"^ 


1000 


0.2 


20.9 


0.11 


D 


2 X 10"" 


1500 


0.4 


13.9 


0.56 


E 


2 X 10"" 


2500 


0.2 


33.5 


0.17 



the right hand side of the equation must reach its max- 
imum value at r=l. Deriving the right hand side of Eq. 
(|14p , one can find the velocity that makes it maximum 
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Since v grows in the inner region, the maximum velocity 
is reached at r = 1, and from Equation 
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7 + 1/ V7+I. 
Continuity (Equation [23]) implies that the value of A is 

1/2 / „, _L 1 \ (37 + 1)7(57+1) 

(26) 
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With the former value of i?, and from Equation 
'7-1 
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the subsonic solution is not attainable (continuity of the 
velocity flow is impossible) and the supersonic branch is 
the only physically viable. In this regime, the flow leaves 
the boundary of the cluster Rc at the local sound speed 
^'max (equal to 1/2 for 7 = 5/3) and is accelerated until 
V = I for r 00. 

Examples of the supersonic flow (velocity and par- 
ticle density) for a group of stars generating different 
values of Massoc , Vw , and Rc , are given in Table 1 . The 
total mass contained up to 10 Rc is also included in the 
Table as an example. A typical configuration of a group 
of tens of stars may generate a wind in expansion with a 
velocity of the order of 1000 km and a mass between 
tenths and a few solar masses within a few pc (tens of 
Rc). We consider hadronic interactions with this matter. 
However, it is to be noted that not all cosmic rays will 
be able to enter the collective wind. The difference be- 
tween an inactive target, as that provided by matter in 
the ISM, and an active or expanding target, as that pro- 
vided by matter in a single or a collective stellar wind, is 
indeed given by modulation effects. The cosmic ray pen- 
etration into the jet outflow depends on the parameter 
e ^ VR/D, where V is velocity of wind, and D is the 



diffusion coefficient, e measures the ratio between the dif- 
fusive and the convective timescale of the particles (e.g.. 

In order to obtain an analytic expression for e for 
a particular star we consider that the diffusion coeffi- 
cient within the wind of a particular star is given by 
([SHHH]) D ^ ^XrC, where is the mean-free-path for 
diffusion in the radial direction (towards the star). The 
use of the Bohm parameterization seems justified, con- 
trary to the solar heliosphere, since we expect that in 
the innermost region of a single stellar wind there are 
many disturbances (relativistic particles, acoustic waves, 
radiatively driven waves, etc.). In the case of a collec- 
tive wind, the collision of individual winds of the par- 
ticular stars forming the association also produce many 
disturbances. The mean- free-path for scattering parallel 
to the magnetic field (B) direction is considered to be 
A|| ~ lOvg — lOE/eB, where rg is the particle gyro- 
radius and E its energy. In the perpendicular direction 
A is shorter, A^ ^ fg- The mean- free-path in the radial 
direction is then given by A^ = A^^ sin^ + A||^ cos^ = 
rg(10cos2 e+siv? e), where cos^'^ e = 1+{B^I B^f .YLeie, 
the geometry of the magnetic field for a single star is rep- 
resented by the magnetic rotator theory ([10]; see also [1] 
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where Vt is the rotational velocity at the surface of the 
star, and B^, the surface magnetic field. Near the star the 
magnetic field is approximately radial, while it becomes 
tangential far from the star, where A^ is dominated by 
diffusion perpendicular to the field lines. This approx- 
imation leads — when the distance to the star is large 
compared with that in which the terminal velocity is 
reached, what happens at a few stellar radii — to val- 
ues of magnetic field and diffusion coefficient normally 
encountered in the ISM. 

Using all previous formulae. 
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Equation (|30[) defines a minimum energy below which 
the particles are converted away from the wind. E"^™ (r) 
is an increasing function of r, the limiting value of the 
previous expression being 

E--ir»R.)^'-^^^(^]^,.S^''* 
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, TeV. (31) 

Therefore, particles that are not convected in the outer 
regions are able to diffuse up to its base. Note that 
i?™'"(r :s> i?*) is a linear function of all i?*, and 

which is typically assumed as ~ 0.1T4o (e.g., PJ)- 
There is a large uncertainty in these parameters, about 
one order of magnitude. The values of the magnetic field 
on the surface of O and B stars is under debate. Despite 
deep searches, only 5 stars were found to be magnetic 
(with sizeable magnetic fields in the range of B^, ^ 100 
G) (e.g., and references therein) typical surface mag- 
netic fields of OB stars are then presumably smaller. 

In the kind of collective wind, we consider that the 
collective wind behaves as that of a single star having a 
radius equal to i?c, and mass-loss rate equal to that of 
the whole association, i.e., Massoc- The wind velocity at 
Rc, K is given by Equation The order of magnitude 
of the surface magnetic field (i.e., the field at i? = Rc) 
is assumed as the value corresponding to the normal de- 
cay of a single star field located within Rc, for which 
a sensitive assumption can be obtained using Equations 
(I28|) and (IMl), 0(10"'^) G. This results, for the whole 
association, in 

'BiRc)\ ( Rc 



0.8 



TeV. (32) 



l^iG ) VO.lpc. 
The value of the magnetic field is close to that typical 
of the ISM, and should be consider as an average (this 
kind of magnetic fields magnitude was also used in mod- 
elling the unindentified HEGRA source in Cygnus, [H]). 
In particular, if a given star is close to Rc its contribu- 
tion to the overall magnetic field near its position will be 
larger, but at the same time, its contribution to the op- 
posite region (distant from it 2 Rc) will be negligible. In 
what follows we consider hadronic processes up to 10 - 
20 i?c, so that a value of the magnetic field typical of ISM 
values is expected. We shall consider two realizations of 
[£;"""(r > 100 GeV and 1 TeV. 



2 7-rays and secondary electrons from a cosmic 
ray spectrum with a low energy cutoff 

The pion produced 7-ray emissivity is obtained from the 
neutral pion emissivity as described in detail in the 
appendix of [15] . For normalization purposes, we use 
the expression of the energy density that is contained 
in cosmic rays, wcr = /g N{E) E dE and compare it 
to the energy contained in cosmic rays in the Earth en- 
vironment, u;cR,e(-^) = -^e(^) ^ '^^i where NpQ is 
the local cosmic ray distribution obtained from the mea- 
sured cosmic ray flux. The Earth- like spectrum, Jq (E) , 
is 2.2E^liJ^ cm-2 s'^ sr'^ GeV'^ (e.g. [ISIIII]), so that 
UJCB.,®{E > IGeV) ~ 1.5eVcm~^. This implicitly defines 
an enhancement factor, as a function of energy 

, , LN(E)EdE , ^ 

^{E) = — ^ ^. . . (33) 



^G-R,®{E) 



We assume that N{E) is a power law of the form N{E) — 
KpE^". Values of enhancement ^ 100 at all energies 
are typical of star forming environments (see, e.g., |15[ 
[TC1|171I191I^U[|^ ) and they would ultimately depend on 
the spectral slope of the cosmic ray spectrum and on the 
power of the accelerator. For a fixed slope, harder than 
that found in the Earth environment, the larger the en- 
ergy, the larger the enhancement, due to the steep decline 
(oc iJ~^-^^) of the local cosmic ray spectrum. In what 
follows, as an example, we consider enhancements of the 
full cosmic ray spectrum (for energies above 1 GeV) of 
1000. With such fixed ^, the normalization of the cos- 
mic ray spectrum, Kp, can be obtained from Equation 
([33]) for all values of the slope. Note that Kp cx and 
thus the flux and 7-ray luminosity, and L^, are lin- 
early proportional to the cosmic ray enhancement. We 
compute secondary particle production (electrons from 
knock-on interactions and electrons and positrons from 
charged pion decay), and solve the loss equation with 
ionization, synchrotron, bremsstrahlung, inverse Comp- 
ton and expansion losses (details are given in the papers 
[201122]). 

We now compute 7-ray fluxes in a concrete example, 
and following Section 1, we consider ^2 Mq of target 
mass being modulated within ^1 pc. The average den- 
sity is ^ 25 cm~^. This amount of mass is typical of 
the configurations studied in Section previously within 
the innermost 20 i?c ~ 2 — 8 pc. To fix numerical val- 
ues, we consider that the group of stars is at a Galactic 
distance of 2 kpc. Using the computations of secondary 
electrons and their distribution, we calculate the 7-ray 
flux when the proton spectrum has a slope of 2.3 and 
2.0. In the latter case, to simplify, we show in Fig.[T]only 
the pion decay contribution which dominates at high en- 
ergies, produced by the whole cosmic ray spectrum. 

The differential photon flux is given by F^{E^) = 
[V/At:D^]Q^{E^) = [M /mp'i-KD'^][Q^{E^)/n], where V 
and D are the volume and distance to the source, and 
M the target mass. In those examples where the vol- 
ume, distance and/or the medium density are such that 
the differential flux and the integral flux obtained from 
it above 100 MeV with the full cosmic ray spectrum is 
greater than instrumental sensitivity, a modulated spec- 
trum with a 100 GeV or a 1 TeV energy threshold might 
not produce a detectable source in this energy range. 
However, the flux will be essentially unaffected at higher 
energy. The left panel of Fig. [T] shows that wind modu- 
lation can imply that a source may be detectable for the 
ground-based Cerenkov telescopes without even being 
close to be detected by instruments in the 100 MeV - 10 
GeV regime (like EGRET or the forthcoming GLAST). 
The right panel of Fig. [1] presents the integral flux of 
7-rays as a function of energy, together with the sensi- 
tivity of ground-based and space-based 7-ray telescopes. 
The sensitivity curves shown are for point-like sources; it 
is expected that extended emission would require about 
a factor of 2 more flux to reach the same level of de- 
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Fig. 1 Diflterential (left) and integral (right) fluxes of 7-rays emitted in a non-modulated and a modulated environment. 
The bump at very low energies in the left panel is produced because we show leptonic emission coming only from secondary 
electrons. Above ~ 70 MeV the emission is dominated by neutral pion decay. Also shown are the EGRET, GLAST, MAGIC 
and HESS sensitivities. Note that a source can be detectable by lACTs and not by GLAST, or viceversa, depending on the 
slope of the cosmic ray spectrum and degree of modulation. Right: Opacities to 77 pair production in the soft photon field of 
an 04V-star at 10, 100 and 1000 i?*, and in the collective photon field of an association with 30 stars distributed uniformly 
over a sphere of 0.5 pc. The closest star to the creation point is assumed to be at 0.16 pc, and the rest are placed following 
the average stellar density as follows: 1 additional star within 0.1, 2 within 0.25, 4 within 0.32, 8 within 0.40 and 14 within 
0.5 pc. 



tectability. Table 3 summarizes these results. From Ta- 
ble 3 and Fig. [1] we see that there are different scenarios 
(possible relevant parameters are distance, enhancement, 
degree of modulation of the cosmic ray spectrum and 
slope) for which sources that shine enough for detection 
in the GLAST domain may not do so in the lACTs en- 
ergy range, and viceversa. 

Finally, In Fig. [1] (rightmost panel) we show the value 
of the photon photon opacity, t{Ej), for different pho- 
ton creation sites distant from a 04V-star 10, 100, and 
1000 i?0, with = 12i?Q and T^s = 47400 K. Unless a 
photon is created hovering the star, well within 1000 i?*, 
7-ray opacities are very low and can be safely neglected. 
This is still true for associations in which the number 
of stars is some tens. Consider for instance a group of 
30 such stars within a region of 0.5 pc (the central core 
of an association). The stellar density is given by Eq. 
([T]); and the number of stars within a circle of radius R 
progresses as A/" = N{R/Rc)^. Fig. [T] shows the collec- 
tive contribution to the opacity in this configuration is 
also very low, since the large majority of the photons are 
produced far from individual stars. However, this is not 
the case if one considers the collective effect of a much 
larger association like the center of Cygnus OB 2 ([23]). 
Reimer demonstrated that even when a subgroup of stars 
like the ones considered here is separated from a super 
cluster like Cygnus 0B2 by about 10 pc, the influence of 
the latter produces an opacity about one order of mag- 
nitude larger than that produced by the local stars. But 
even in this case, Fig. [1] shows that this opacity is not 
enough to preclude escape from the region of the local 
enhancement of stellar density. 



3 Conclusions 

We have studied collective wind configurations produced 
by a number of massive stars, and obtained densities and 
expansion velocities of the stellar wind gas that is target 
for hadronic interactions in several examples. We have 
computed secondary particle production, electrons and 
positrons from charged pion decay, electrons from knock- 
on interactions, and solve the appropriate diffusion-loss 
equation with ionization, synchrotron, bremsstrahlung, 
inverse Compton and expansion losses to obtain expected 
7-ray emission from these regions, including in an ap- 
proximate way the effect of cosmic ray modulation. Ex- 
amples where different stellar configurations can produce 
sources for GLAST and the MAGIC/HESS/ VERITAS 
telescopes in non-uniform ways, i.e., with or without the 
corresponding counterparts were shown. Cygnus OB 2 
and Westerlund 1 maybe two associations where this sce- 
nario could be at work ([22]). 
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